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We derive the effective action of the hard large- a; valence charges up to fourth order in their density. 
Such non-Gaussian weight functionals contribute at leading order in Nc to the connected two-gluon 
production diagrams which determine di-hadron correlations. The corresponding diagrams are not 
necessarily (highly) suppressed by the density of valence charges since their infrared divergences 
differ from those obtained in a Gaussian theory. Therefore, it appears prudent to include such 
higher dimensional operators when determining initial ensembles for JIMWLK evolution of higher 
^ — ^ . n-point functions of Wilson lines. 

(N 

>>■ I. INTRODUCTION 

(Multi-) Particle production cross sections in hadronic collisions at high energies are related to expectation values 
, of traces of products of Wilson lines which resum multiple scattering and high gluon density effects. The evolution of 
P\| ' such operator expectation values is described by the JIMWLK functional renormalization group equation [J]. In the 
large- A^c limit, and in a Gaussian approximation for the effective action, they reduce to the BK equation j^J for the 
dipole forward scattering amplitude. 

The high-energy evolution equations require initial conditions at a rapidity Y = log xq/x = Q (xq is often assumed 
I . to be about 10"^). They have been derived by McLerran and Venugopalan Q in the limit of an infinitely large 
nucleus. In the MV model, the large-x valence charges act as recoilless sources for the soft, small-x gluon fields. As 
^ , A — ^ oo then, the variance of the "valence" charge density ~ A^^^ grows large and the distribution of color charges 
should be Gaussian, 
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It should be noted that high multiplicity proton-proton collisions may also correspond to unusually high valence 
charge densities (see, for example, ref. [iJ]) and so would also be described effectively by this action with A^s ^ 1- 

Nevertheless, in reality the mass number A resp. the number of valence charges is finite, in particular in (high- 
multiplicity) pp as well as peripheral AA collisions. It is therefore interesting to consider extensions of the MV-model 
action involving higher powers of the color charge density. In fact, Jeon and Venugopalan have derived [5| an 
"odderon" operator contribution —d'^'^'^p'^p^p'^/Ks, where the cubic coupling K3 ^ g'^{A/nR^)^ . Below, we shall show 
that at quartic order in p a contribution 6°'^ 5'^'^ p"" p^ p'^ p'^ / K4, to the effective action arises, with K4 ~ g'^{A/ -kR^Y . 

_ ^ ^ Our motivation derives from the recent observation of a "ridge" in two-particle correlations from high-multiplicity 

] pp as well as heavy- ion collisions [7|- The ridge refers to a correlation which extends over several units^ of (relative) 
$^ : rapidity for particles with similar azimuthal angle, A0 ^ tt. For high-multiplicity pp collisions at ^/s — 7 TeV the 
CMS collaboration found that the anxplitude of the "ridge" correlation peaks about transverse momenta on the order 
of 2-3 GeV, a semi-hard scale. Ref. Q argued that within the framework of high energy evolution this effect would 
correspond to production of two small - a; g luons with relative rapidity Ay > 1 from two evolution ladders which 
connect to the same large- a; sources [13, 

With a Gaussian action one finds that the two-particle distribution minus the product of two single-particle distri- 
butions, 
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It is presently not known whether the correlation simply satisfies boost invariance or whether it diminishes for rapidity intervals on the 
order of l/a^ [EI3|- 
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is suppressed by a factor 1/{N^ — 1) relative to the uncorrelated part j^'fj^y y \ d^qdy /• subleading order 

in A^c, however, Gaussian factorization of the four-point function which determines ^ is violated by JIMWLK 
evolution ^121] . 

This problem could be cured by a dynamically generated correlation length (in the transverse plane) of small-x 
gluon fields [l^ . For a discussion of how color charge correlations develop as a hadron or nucleus is boosted to high 
rapidity we refer to ref. [TJ] . Averaging of n-point functions of Wilson lines with a local Gaussian would not apply 
when those n points are within one correlation length of each other and consequently a connected contribution should 
arise at leading order in Nc [l^ ■ This can be seen also from a model which assumes splitting of the four-point function 
into two-point functions ("ladders") with non-zero individual transverse momenta [15 1. 

Here, we show that a non-Gaussian initial distribution of valence charge naturally arises when A < oo and that it 
leads to particle correlations at leading order in Nc- Furthermore, that the infrared behavior of the leading connected 
two-particle production diagram is different from the case of a quadratic action. This fact may affect the relative 
A-dependence of the correlations resulting from the quadratic vs. the higher-order terms in the action. 

We focus on the kinematic regime where the transverse momenta p, q of the two produced gluons are somewhat 
larger than hut on the order o/the scale where evolution in rapidity is non-linear. For a proton, this so-called saturation 
scale Qs{x) is expected to be on the order of 1 GeV at x ~ a few times lO"**, on average; high-multiplicity collisions 
should correspond to configurations with significantly higher parton densities. For nuclei, the valence charge density 
should be boosted by a factor of ^ A^^^. In this kinematic regime, the transverse momenta of the produced gluons 
are, to a significant part, due to the intrinsic transverse momentum from the small-a: evolution ladders. Further, the 
relative azimuthal angle of p and q is taken to be small, ^ tt. 

When logp^/Qs 3> 1 and logg^/Q? 3> 1 the situation is different. This case has been analyzed within the double- 
logarithmic approximation in ref. [la]. They show that here a single "flip" (or "recombination") between the two 
evolution ladders is suppressed and that two such flips are required, one below and the other above the rapidities of 
the produced dz-jets. At very high energies such flips between ladders may occur at a hard transverse momentum 
scale not too far from the hard dijet vertices themselves [T6| . 



II. DERIVATION OF THE EFFECTIVE ACTION BEYOND QUADRATIC ORDER 

In this section we derive the form of the effective action for a system of k quarks in SU(3) following the methods 
developed in [5|. The probability to find this system in a representation labeled by the integers (m, n), also denoted 

{k) 

as Dynkin labels, is given by dmn Nrn'n where 
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We are interested in the representation with the largest weight and so consider the limit fc ^ r7i,n 3> 1. Thus, we 
expand the factorials in each G in powers of 1/fc up to and including order ^ -p.^ xhe leading terms of the form 



t — - can be written as 
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where C2, C3, C4 are the SU(3) Casimir operators for the representation (m, n), given by (see appendix) 
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Before we rewrite the action S{m, n; k) in terms of the color charges p, it is perhaps useful to remind the reader of 
this derivation and to illustrate it in the simpler case of SU(2) spins [5|]. In this case, the representation R is labeled 
by one index, I, and there is only one independent Casimir defined as C2(i?)l/? = S"-'' LJj = Lj^, where are the 
generators of the representation R. For a large representation we have C2 — ({£ + 1) ~ £^ as the quadratic Casimir. 
One can then write £^ = £i £i where £i is a vector describing spin and express the action in terms of invariants formed 
by multiplying £i's (in this example, there is only one invariant, £'^) as S ~ Nc£'^ /k. For large £, one introduces a 
classical charge density per unit transverse area, = g£''/ A'^x. Next, k is expressed in terms of number of valence 

quarks in a nucleus as k = NcAA'^x/nR'^ which then leads to 



C,,{R) = —F^"- tvT^' ■■ - T^- . (11) 



with /i2 EE 

We start by defining the Casimirs for general Nc in terms of the generators T^' of the representation R, 

C„(i?)llfl = F^^' - T^' • • • T^" , (10) 

where 1^ is the unit matrix in the representation R. Taking the trace of both sides gives Cn{R), the n*"^ Casimir of 
the representation R: 

1 

The color tensors F^^' " are the most general color invariant tensors one can construct out of the SU(iVc) orthonor- 
mal basis tensors. For example, F'^'"'^ = for n = 2, while for n = 4 these basis tensors are given by 5°"^ 5'^'^ plus 

permutations, c?"*"^ (T'^'^ plus permutations, and d"'"^ f'^'^'^ plus permutations. More explicitly, 

C2{R) = ^^^'HrTj^T^ (12) 

Since the number of independent Casimirs of S\J{Nc) is equal to its rank, A'^ — 1, there are only 2 independent Casimirs 
for SU(3). These are usually taken to be C2 and C3. This can be seen explicitly by noting that for Nc ~ 3 one has 
the constraint 

^abe^cde ^ ^ace^bde ^ ^ade^bce ^ }_ (^gabgcd _^ gacgbd _^ gadgbc^ (^^g^) 

3 

and 3a + /3 = 1 so that the quartic Casimir can be written as the square of the quadratic invariant. To proceed, we 
write the action in terms of a — 1 dimensional vector related to the second Casimir via |Q| = y/Q'^Q"' = VC2', 
the vector is analogous to the angular momentum vector in our SU(2) example from above: 



S: 



We then follow McLerran-Venugopalan and define the color charge per unit area = gQ°'/A^x to finally arrive at 

"(J'^Vp^ d"'"' P" S'^^Scd ^gac^bd ^^ad^bc ^ ^ ^ ^ 



S[p{x)] ~ / d^x 



-p p p p 



(14) 



K3 K4 

To write the action in this form we have used eq. and k — NcA^^. This assumes that the large-x sources 
correspond to NcA valence quarks. On the other hand, if initial conditions for small- a; evolution are set at, say, 
xo = 0.01 then the number of "valence" charges would be bigger although still parametrically proportional to A and 
to Nc. 

The couplings in this action are given by, parametrically, 

^ 0(3^^1/3) , (15) 
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0{g'A'/^) , (16) 



«4 ^ 18 (f]^ - 0(/A) . (17) 
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The non-Gaussian terms in the action thus involve additional inverse powers of gA^^^. Also, the expressions from 
above are obtained by averaging over the transverse (impact parameter) plane; at a fixed distance b from the center 
of a nucleus one should replace A/ttR'^ by the thickness function T{b). 

The action (|14p is coupled to the soft gauge fields f]\ . At leading order and in the classical approximation this leads 
to the relation ([22)1 below. 

III. TWO-POINT FUNCTION AND UNINTEGRATED GLUON DISTRIBUTION 




FIG. 1: Two-point function at leading order in I/K4. 



In order to compute the contribution of the quartic term to a physical observable, such as single-inclusive hadron 
production in the forward region (Tt} , one needs to perform the color averaging of the dipole operator with this new 
action: within the Glauber- Mueller approach the expectation value of the full dipole operator (tr VxV^), where V 
denotes a Wilson line, resums multiple scattering effects on a dense target 0. This can in principle be done numerically 
using lattice gauge theory methods [18]. Nevertheless, it is useful to consider a limit where this expectation value can 
be evaluated analytically, such as the dilute limit of the dipole cross section given by the two-point function of color 
charges. Even then, one can not perform the integration over p analytically when the action is not quadratic. We 
therefore resort to a perturbative expansion in I/K4 and keep only the first term in the expansion of the exponential 
of the quartic term. We then compute the two-point function of the color charge density of hard sources which is 
related to the unintegrated gluon distribution. To first non-trivial order in I/K4, 
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Here, So denotes the Gaussian action. For the two-point function O — p'^Py, the possible contractions are shown 
diagrammatically in fig. [T] for the denominator of eq. (jl8p one amputates the points x and y. We compute the 
functional integral in lattice regularization, i.e. we approximate the two-dimensional transverse space by a lattice with 
Ns sites of area A^x. The two-point function becomes^ 
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K4 ^ c 
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(19) 



The disconnected contribution exhibits both a UV (A^a; — >■ 0) as well as a IR {Ng — > 00) divergence but appears both 
in the numerator and in the denominator and so cancels as usual. The third term in the numerator is due to the 
tadpole diagram from fig. [T] which renormalizes the "bare" parameter p^: 



p' 



1 



sP 



1 



K4 



(20) 



This rcnormalization of p?' will be applied to all diagrams to absorb the insertion of a tadpole into any line. 



The right-hand-side of this expression is to be understood in lattice notation: x and y are discrete points and S^y is a Kronecker symbol. 
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Thus, the two-point function now reads (in continuum notation) 

iplpD^fi'S^'Six-y) o {p*'^{k)p\k'))=p'S'^'{27T)H{k~k') . (21) 

To relate (PT|) to the unintegrated ghion distribution we note that at leading order and in covariant gauge the field 
generated by a fast particle moving in the positive z-direction is related to the charge density by 

A^(x) = (5^+a(a;) = g S''+6ix-)^p{x) , (22) 

where again x denotes a transverse coordinate. This field also satisfies = and so the only non-vanishing field- 
strength is F'^'^ — —d^a. In momentum space we have the relation k^a{k) — gp{k). Thus, we define the unintegrated 
gluon distribution $(fc^) ^ fc^ via 

(p-(fc)p^(fc')) - --p—{27:rS{k - k') $(fc2) . (23) 
^ ' as — 1 

With this convention, ^(fc^) = as{N'l - l)M^/(27r). 



IV. FOUR-POINT FUNCTION 
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FIG. 2: Four-point function at leading order in 1/ hn. 
The four-point function to first non-trivial order in I/K4 is given by 



The different type of contractions which arise at this order are shown in Fig. [21 Eq. ([M)) then gives 

1 
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(25) 



The second line in this expression corresponds to the sum of the first three diagrams shown in Fig.[2J the second term 
due to the disconnected diagram again cancels against the overall normalization factor once the latter is expanded to 
leading order in I/K4. Also, the third term from the second line is again absorbed into the rcnormalized value of p^ 

as in (f20| . In fact such a factor of 1 — 16^ (iV^ -|- l)-^!^ appears at order I/K4 from diagrams of the type shown in 
Fig. |3] which induce the shift p^ — )■ p^ in the last line of eq. (p5)) . Thus, the four-point function to order I/K4 finally 
becomes (switching again from lattice to continuum notation) 



/ a b c d\ ~4 
{PxPvPuPv) = M 
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FIG. 3; Diagram at order I/K4 which renormahzes /i^ in the connected contribution to the four-point function. 



In momentum space, 

{2^tp' 
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(27) 



V. CORRELATED TWO-GLUON PRODUCTION 
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FIG. 4: Correlated production of two gluons with small relative azimuth, <^ tt. 



We now turn to the color factor for correlated production of two gluons as illustrated in Fig. |31 We focus, in 
particular, on the so-called "ridge" kinematics where the gluons are separated by a rapidity interval A77 > 1 but 
where their transverse momenta are nearly parallel so that their relative azimuth <^ tt. We assume that in this 
kinematic regime the dominant contribution arises from a process where the small- a; gluons are produced from two 
separate ladders which connect to the same hard sources, though fToj. 

It was pointed out in ref. ^12i] that the cross section for this process involves a product of two four-point functions 
(projectile and target, respectively) such as the one discussed above. Further, that for a Gaussian action the correlated 
contribution which is left over after one subtracts the square of the single inclusive cross section is suppressed by 

The diagram from Fig. |3] corresponds to [12| 



/eaa'/e'bb'/ecc'/eMd'(p*"(fc2)p*'(fc4)p"(fcl)p''(fc3))(p*"'(p-fc2)p*''(g- W'b-^O/l'Z-fca)) 



(28) 



with e and e', respectively, the color indices of the two produced gluons; the ladder momenta ki will be integrated 
over. The (squared) single-inclusive cross section is obtained from 



feaa'fe'bb'fecc'fe'dd'{p*''{k2)p%k,))(^p*'{k4)p''{k3)) (p*''\p- k2)p'\p^ k^)) ^ 
{2TTf Nl{Nl - if (TTR^f p"" d{ki - k2) S{k3 " fc4) . (29) 
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Translational invariance in the transverse plane brings along a factor of S{k = 0) = 7ri?^/(27r)^ for each of the two 
produced gluons. Using ([15]) the prefactor is of order 

(single inclusive)' ^ N',{N^ - 1)' {^R^f (f]^)' ^ (.9^'/')' (^^c " 1)' ■ (30) 




FIG. 5: Two-gluon production diagram which gives an angular coUimation about A<ji = over several units in relative rapidity, 
|j/p ~ J/?! ^ 1 (the so-called "ridge") [^. Such diagrams arise from a Gaussian action and only involve the two-point function 
(unintegrated gluon distribution); they are suppressed by ~ 1/{N^ — 1) relative to uncorrelated production. 



The remaining contractions in eq. (|28p generate connected diagrams corresponding to correlated production. In the 
limit K4 — >■ 00 they are suppressed by a factor 1/{N^ ~ !)■ For example, for the diagram from fig. O 

feaa'fe'wfecc'fe'dd'{p*^{k2)p\h))(^p*\h)p-{h)) (p* (p - fc^ )/ (p " fcl ) ) {p* {Q ~ h) p''' {Q ~ h)) 

= {2TTf Nl{Nl - 1) TiR^ fl^ 6{ki - fca) 6{ks - h) 5{ki - k^) . (31) 



Due to the additional ^-function as compared to eq. ((29|) there is, of course, only one single independent ladder 
momentum to integrate over. This diagram will therefore depend on the relative azimuthal angle between p and q; 
it is proportional to 



Fig. 5 ^ NliNl - 1) .R^ ^ j . (32) 

Its contribution is therefore maximal at — and minimal at A(/) = 7r/2 [8]. Also, the integral is quadratically 
divergent from the region fc^ ^ p',q' and k ~ p ~ q. If a fixed non-perturbative cutoff is used then this 
contribution \s ^ , i.e., suppressed by ~ A~'^/^ as compared to uncorrelated production, eq. ([50]) . On the other 
hand, if such cutoff is in fact provided by the non- linear saturation scale [l3|, Qs ^ A^/^ , then (|31l) is of order ^ A^/^ , 
down by A~^ as compared to uncorrelated production'^. 

For finite K4, however, we can see from eq. (I26I27P that the term ^ 5°''^ 5^'^ produces the same overall color factor as 
the square of the single-inclusive cross section. The leading order (in I/K4) connected contribution of order N^{N^ — 1)'^ 
to (|28p is shown diagrammatically in fig. [6] and is given by'^ 

-2(27r)« A^/i4 S''^5'^''S'''^'S''''''f,aa'fe'bb'fecc'fe'dd' S{ki - ^2) 6{k3 - ki) 5{ki + A^g - ^ ^^4) 
TT^ K4 

= -(27r)6 4 Nc{Nl - If ^ 5{ki ~ k2) S{k3 - ki) . (33) 



dp 1 1 



^ Note that in practice this ~ 1 /A suppression of the correlated contribution is usually removed by multiplying with the multiplicity per 
unit rapidity, dN/dy. 

The overall factor of 2 accounts for projectile target flip. 
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FIG. 6: Diagram for correlated two-gluon production from the quartic action. This contribution is of the same order in g and 
A as that from fig. [5] but enhanced by a factor of A'^^ — 1. 



The contribution from this diagram to the two-particle spectrum is 
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(34) 



Thus, this contribution is independent of the azimuthal angle A(f) between p and q. Also, its sensitivity to the infrared 
cutoff is only logarithmic and would therefore not modify the ^-dependence of the prefactor. 
From then the prefactor of ([55]) is of order 



Fig. 6 ^ N^{N! 



1 firR' 



A 



^8^5/3 ^2(^2 „ 



(35) 



This is of the same order in and Nc as the uncorrelated contribution ([5D|) but suppressed by one power of the mass 
number A. However, it is more relevant to compare the contribution from the quartic interaction to the correlation, 
eq. (|35l) . to that from the Gaussian action, eq. (|32l) . If the latter is cut off at a constant, A- independent scale 
then it is parametrically enhanced by a factor of A-^^^ and the corrections from the quartic term in the action can 
be computed perturbatively (near A = oo) by an expansion in I/K4 (as we have done). One may expect that such 
corrections are large when A^^^ ~ 6 since they are enhanced by a factor of — 1. 

On the other hand, if the infrared cutoff ~ ~ A^/^ then the non-Gaussian action generates "corrections" to 
two-particle correlations which appear at the same order in A as the contribution from the Gaussian action (up to 
logarithms of A); in addition, those "corrections" are enhanced by a large color factor — 1. This case requires a 
non-perturbative calculation of the correlation to all orders in I/K4. 





FIG. 7: Connected contribution to the four-point function at order 



We conclude this section by noting that the "odderon" operator can also contribute to the connected two-particle 
production cross section at order ^ 1/k§. This contribution is obtained by expanding the exponential of the odderon 
operator to second order. Performing the averaging of the four-point function shown in fig. [7] over the hard color 
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sources like in eq. (|24|) results in five Wick contractions and is of order / k"^. Multiplying by the leading ^ Kg 
contribution from the target side gives, in all: 

+ nR' Nl{Nl -l)^N.-±-yt^ ^ + nI{NI - 1) [n. - ±) . (36) 

This is of the same order in A as the contribution ~ I/K4 from eq. psp although it is suppressed by one power of N^- 
It appears sensible that numerical solutions of JIMWLK evolution also include the odderon operator in the action 
for the initial condition. On the other hand, terms of order ~ do not matter for the four-point function (they do 
matter for the six-point function and higher), they simply renormalize K4. 



VI. SUMMARY 



In summary, we have computed corrections to the quadratic MV-model action up to quartic order in the density of 
hard (large-x) sources. Such terms are accompanied by additional powers of 1/ gA^/^ . This action could be employed 
to determine initial conditions for the JIMWLK high-energy evolution of various n-point operator expectation values. 
We have motivated this by showing that at leading order in the quartic coupling this action gives a connected 
contribution to two-gluon production (i.e., to correlations) already at leading order in TVc, enhanced by a factor of 
— 1 over the connected contribution from a (local) quadratic action. Such non-Gaussian actions may also be 
important for applications to dijet production with transverse momenta not very far above the saturation scale [l9| . 



Acknowledgments 

We thank S. Jeon for extensive discussions on the derivation of the MV action in ref. 's'l and A. Kovner for 
communications which initiated this work. We gratefully acknowledge support by the DOE Office of Nuclear Physics 
through Grant No. DE-FG02-09ER41620, from the "Lab Directed Research and Development" grant LORD 10-043 
(Brookhaven National Laboratory), and from The City University of New York through the PSC-CUNY Research 
Award Program, grants 63382-0041 (A.D.) and 63404-0041 (J.J-M.). 



VII. APPENDIX 



In this appendix we illustrate how one can express the values of various Casimirs in terms of the Dynkin labels n 
and m of SU(3). The representation {171,11) is reached by multiplying the (1,0) fundamental representation k times. 
The corresponding Young tableau thus has k boxes in = 3 rows. The first row has hi =m + n+{k — m — 2n)/3 
boxes; the second row has h2 = n + {k — m — 2n)/3 boxes; and the last row has = {k ~ m — 2n)/3 boxes. 

It is useful to define the set n- for the U(N) group from these "row lengths" hi via n- — hi + {N + l)/2 ~ i: 

, k — m — 2n k + 2m + n , ^ 

n[ = m + n+ + 2 - 1 = + 1 (37) 

, k — m — 2n k — m ^ n , ^ 

n'^ = n+ + 2 - 2 = (38) 

,fc — 771— 2n k — m — 2n 
«3 = 5 + 2 - 3 = 1 (39) 



n[ = k . (40) 



To pass from U(3) to SU(3) one needs to shift rii = n'i — J2j n'j/N to get 

2m + n , ^ 

ni = ^— + 1 (41) 

— m + n , , 

n2 - (42) 

—TO — 2n ^ , 

713 = 5 1 • (43) 
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Therefore, 

y^n^ = = - (m^ + + mn) + 2{m + n + 1) (44) 

O 

I 

= - + +71^ + 3n + 3m + mn)^ . (45) 



The quadratic casimir is constructed from Ui via 

12 3 



C2(m-, n) = ^ — = - {m^ + 11? + mn) + (m + n) . (46) 



For a single quark corresponding to the fundamental representation (m = 1, n = 0) of QCD, we get C2 = {N^ — 
l)/(2iV,) = 4/3. 

To compute the cubic Casimir we require the set of h — nt + (N — l)/2, 
2m + n ^ 

k = — + 2 (47) 

h = ^ + 1 (48) 

—m — 2n , 
^3 = ^ (49) 

^ N (50) 

i 

2 

^Y^il = iV + ^71^ = - (m^ + 71^ + m7i) + 2(m + n) + 5 (51) 

i i 

E, , , , , 2771 + n Am? + r? -\- Amn 3m + 6n — nm + — m? 
hlj = h{l2 + h) + l2h = 2 (52) 



3 

771^ + r? + mn + 3m + 3n 



(53) 

E,o 277?,'^ — 277'^ + 'im^n — 3?7777^ TO . ^ 

= h 3771"^ + 77^ + 277777 + 7777 + 577 + 9 . (54) 

i 

These can then be used to construct the cubic Casimir following [23| and lead to 

(73(771, n) = ^ (777 + 277 + 3) (77 + 2777 + 3) (m-n) . (55) 
18 

Constructing C4{m,n) requires the vectors 77^ and can be done similarly to (72(777,77); this leads to 

1 2 

(74(777, 77) EE - (777'' + 77"* + 277777^ + 27^7^77 + 3777^77^) + - (m^ + 7t^ + 2777^71 + 277777^) 

y o 

+ i (5771^ + 577^ + ll7777j) - i (tR + 77) . (56) 
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